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We analyze the effect of local decoherence of two qubits on their entanglement and the Bell in- 
equality violation. Decoherence is described by Kraus operators, which take into account dephasing 
and energy relaxation at an arbitrary temperature. We show that in the experiments with super- 
conducting phase qubits the survival time for entanglement should be much longer than for the Bell 
inequality violation. 
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Entanglement of separated systems is a genuine quan- 
tum effect and an essential resource in quantum informa- 
tion processing^ Experimentally, a convincing evidence 
of a two-qubit entanglement is a violation of the Bell 
inequality^ in its Clauser-Horne-Shimony-Holt^ (CHSH) 
form. However, only for pure states the entanglement 
always- results in a violation of the Bell inequality. In 
contrast, some mixed entangled two-qubit states (as we 
will see, most of them) do not violate the Bell inequality, 5 
though they may still exhibit nonlocality in other ways. 6 
Distinction between entanglement and Bell-inequality vi- 
olation, in its relevance to experiments with supercon- 
ducting phase qubitSfi is the subject of our paper. 

The two-qubit entanglement is usually characterized 
by the concurrence 8 C or by the entanglement of 
formation^ 9 - which is a monotonous function 8 - of C . Non- 
entangled states have C = 0, while C = 1 corresponds to 
maximally entangled states. There is a straightforward 
way 8 - to calculate C for any two-qubit density matrix p. 
The Bell inequality in the CHSH form 3 is |5| < 2, where 
S = E(a, b) - E(a, b') + E(a', b) + E(af, b') and E(a, b) 
is the correlator of results (±1) for measurement of two 
qubits (pseudospins) along directions a and b. This in- 
equality should be satisfied by any local hidden-variable 
theory, while in quantum mechanics it is violated up 
to | if?! = 2^/2 for maximally entangled (e.g., spin-zero) 
states. Mixed states produce smaller violation (if any), 
and there is a straightforward wayi^ to calculate the max- 
imum value S+ of IS* | for any two-qubit density matrix. 

For states with a given concurrence C, there is an exact 
boundii for S+: 2\/2C < S+ < 2vT + C 2 (we consider 
only S+ > 2), so that the Bell inequality violation is guar- 
anteed if C > l/y/2. For any pure state the upper bound 
is reached: S+ = 2yl + C 2 , so that non-zero entangle- 
ment always leads to S+ > 2. The distinction between 
entanglement and Bell inequality violation has been well 
studied for so-called Werner states 5 - which have the form 
p = fp s + (1 — /)p m ix, where p s denotes the maximally 
entangled (singlet) state, and p m i X = 1/4 is the density 
matrix of the completely mixed state. The Werner state 
is entangled for— / > 1/3, while it violates the Bell in- 
equality only wheni^ / > 1/ y/2. 

The Werner states, however, are not relevant to 
most of experiments (including those with supercon- 



ducting phase qubits^), in which an initially pure 
state becomes mixed due to decoherence (Werner states 
are produced due to so-called depolarizing channel^). 
Recently a number of authors have analyzed ef- 
fects of qubit decoherence on the Bell inequality 
violation 1 ^ 1 ^! 5 .^ 6 and entanglement ^^^^23 
Best-studied models of decoherence in this context are 
pure dephasin g 12 ' 13 ! 15 ' 19 ! 21 ' 23 and zero-temperature en- 
ergy relaxation ; 14 ' 16 ' 18 ' 22 while there are also papers 
considering a combination of these mechanisms ) 17 ' 20 
high-temperature energy relaxation? 1 - 4 , and non-local 
decoherence . 14 ' 15 ' 23 In particular, for the case of pure 
dephasing it has been show n 19 ' 20 that the concurrence 
C decays as a product of decoherence factors for the 
two qubits, and therefore a state remains entangled for 
arbitrarily long time; moreover, the calculation of S+ 
show o 12 ' 13 that the Bell inequality is always violated also. 
For the case of zero-temperature energy relaxation it has 
been shown that entanglement can still last foreve r 16 ' 18 ' 22 
(depending on the initial state), while a finite survival 
time has been obtained^ 6 - for the Bell inequality viola- 
tion. 

In this paper we consider a two-qubit state decoher- 
ence due to general (Markovian) local decoherence of 
each qubit (including dephasing and energy relaxation 
at a finite temperature) and assume absence of any other 
evolution. For this model we compare for how long an 
initial state remains entangled (C > 0), and for how long 
it can violate the Bell inequality (S + > 2). In particular, 
we show that for typical (best) present-day parameters 
for phase qubits 7 - these durations differ by ~ 8 times. 

Before analyzing this problem let us discuss which frac- 
tion of the entangled two-qubit states violate the Bell 
inequality. This question is well-posed only if we intro- 
duce a particular metric (distance) and corresponding 
measure (volume) in the 15-dimensional space of den- 
sity matrices. Various metrics are possible; let us choose 
the Hilbert-Schmidt metric? 1 -^ for which the geometry 
in the space of states is Euclidean. Then random states 
p with the uniform probability distribution can be gen- 
erated as- 24 p = A/tr(A <! A), where A is a 4 x 4 matrix, 
all elements of which are independent Gaussian complex 
variables with the same variance and zero mean. Us- 
ing this method, we performed Monte-Carlo simulation, 
generating 10 9 random states and checking if they are 
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entangle d 25 ' 26 ! 27 and if they violate the Bell inequality. 10 
In this way we confirmed that 75.76% of all states are 
entangled 2 - and found that only 0.822% of all states vio- 
late the Bell inequality. Therefore, only a small fraction, 
1.085% of entangled states violate the Bell inequality. 

Now let us discuss the effect of decoherence. For one 
qubit it can be described by the Bloch equations 2 ^ (we 
use the basis of the ground state |0) and excited state 
|1)) and characterized by the energy relaxation time Ti, 
dephasing time T 2 (T 2 < 2T\ ) , and the Boltzmann factor 
h = exp(— A/#), where A is the energy separation of the 
states and 8 is the temperature. The usual solution of the 
Bloch equations can be translated into the language of 
time-dependent superoperator C for the one-qubit den- 
sity matrix p, so that p(t) = £[p(0)] = Ej=i KiPityKj > 
where four Kraus operators Ki can be chosen as 



K, = 








X/VT~g 



K 4 = 



Vhg , (1) 
' ' 1 > 



where 9 = [1 - exp(-t/Tj.)]/(l + h), A = exp(-t/T 2 ), 
and in our notation |1) = (1,0) T , |0) = (0, 1) T . It is 
easy to check that the term under the square root in 
K3 is always non-negative and equals (for t > 0) only 
if T 2 = 2T\ and 6 — 0. Notice that choice of the Kraus 
operators Ki is not unique (though limited to the unitary 
freedom of quantum operations 1 -) and, for instance, the 
Kraus operators presented in Ref. [l] for the special cases 
of depolarizing channel (Ti = T%, = 00) and energy 
relaxation (T 2 = 2T X ) differ from Eq. (JIJ. 

In general, decoherence of two qubits is described 
by many parameters (out of 240 parameters describing 
a general quantum operation only 15 parameters de- 
scribe unitary evolution). We choose a relatively sim- 
ple but physically relevant model when the decoher- 
ence is dominated by local decoherence of each qubit. 
(Non-local decoherence would be physically impossible 
in the case of large distance between the qubits.) The 
model now involves six parameters: T"' 6 , , and 
h ab — exp(—A at b/6 a> b), where subscripts (or super- 
scripts) a and b denote qubits, and the evolution is de- 
scribed by the tensor- product superoperator C — C a ® C b 
(which is completely positive because of complete posi- 
tivity of C a ,b)- This superoperator contains 16 terms: 

P (t) = c[p(Q)} = Et,-=i*«p(°)*£-> K a = Kf^K), 

where operators K^' b are given by Eq. {T]) for each qubit. 
As an initial state we consider an "odd" pure state 
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cos /3 1 10) + e ia sin/3|01) 



(2) 



(0 < j3 < 7r/2), which is relevant for experiments with 
the phase qubits.— Since the parameter a corresponds to 
z-rotation of one of the qubits, while decoherence as well 
as values of C and S+ are insensitive to such rotation, all 
results of our model have either trivial or no dependence 



on a. The evolution of the state ([2]) due to local decoher- 
ence C can be calculated analytically, and at time t the 
non-vanishing elements of the two-qubit density matrix 
p are 

pu(t) = (1 - g a )h b g b cos 2 f3 + h a g a (l - g b ) sin 2 f3, 
P22{t) = (1 - 5a)(l - h b g b ) cos 2 13 + h a g a g b sm 2 (3, 
P33(t) = g a h b g b cos 2 (3 + (1 - h a g a )(l - g b )sm 2 (3, 
Pa(t) = g a (l- h b g b ) cos 2 (3 + (1 - h a g a )g b sin 2 0, 
P32(t) - p* 23 (t) - exp(-i/T 2 a - i/7 1 2 b )e M (sin2/3)/2, (3) 

where g a b are defined below Eq. (fTJ), and p t j sub- 
scripts i,j = 1,2,3,4 correspond to the basis 
{|11), 1 10), 1 01} , 1 00)}. These equations become very sim- 
ple at zero temperature because then h a = h b = 0. No- 
tice that the dephasing times T^'" enter Eqs. §5§ only 
through the combination 1 /T 2 + 1 /T| (this is not so for a 
general initial state) , so that the two-qubit dephasing can 
be characterized by one parameter T 2 = 2/(l/T 2 +1/T|). 
For the state © the concurrence i o 14 ' 20 

C = 2max{0, \pm\ - Vpilpil}, (4) 
and the Bell inequality parameter S+ i o 10 ' 16 

S+ = 2max{2V2|p 23 |, ^Itol 2 + (1 - 2 Pll - 2p 4 4) 2 }, 

(5) 

while for the initial state C = sin 2/3 > and S + — 
2Vl + C 2 > 2. Notice that the first and second terms 
in Eq. ([5]) correspond to the "horizontal" and "verti- 
cal" measurement configurations, using the terminology 
of Ref. H3- Equations ((3]), ((4j), and (0 are all we need to 
analyze entanglement and Bell inequality violation. 

Notice that for a pure dephasing (T" = T\ = 00) we 
have pn = P44 = 0, and therefore 



C = cxp(-2i/T 2 ) sin 2/3, S+ = 2y/l + C 2 . (6) 

In this case at any t the state remains entangle d 19 ' 20 and 
violates the Bell inequalit y 12 ' 13 (It also remains within 
the class of states producing maximal Bell inequality vi- 
olation for a given concurrence^ 1 -) In the case when both 
dephasing and energy relaxation are present but temper- 
ature is zero, 8 a = 9 b = 0, the concurrence C is still 
given by Eq. © and lasts forever j 16 ' 22 however S + docs 
not satisfy Eq. © and, most importantly, the Bell in- 
equality is no longer violated after a finite timeJS- Fi- 
nally, in presence of energy relaxation at non-zero tem- 
perature (at least for one qubit) the entanglement also 
vanishes after a finite time, as seen from Eq. in which 
lim^oo pnp44 7^ 0. 

Let us consider in more detail the case when both 
dephasing and energy relaxation are present, but tem- 
perature is zero and Tf = T x b = T x . Then Eq. © 
for S+ becomes very simple since pn = and P44 = 
1 — exp(— t/Ti). The time dependence S+(t) consists of 
three regions: at small t it is always determined by the 
second term 3 ^ in Eq. (©, then after some time t\ the first 
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FIG. 1: The two-qubit entanglement duration te in units 
of the dephasing time T2 for the maximally entangled ini- 
tial state (P — 7r/4) and several values of the temperature 8. 
Dashes lines correspond to Eq. (|8]). 



FIG. 2: The duration tb of the Bell inequality violation (as- 
suming p = tt/4) for A/6 = 15 (solid line) and A/0 = 
(dotted line). The dashed line: t b /T 2 = ln[Ti/(4r s )]/4. 



term becomes dominating, while after a later time t 2 the 
second term becomes dominating again. Notice that in 
the second region S+ = 4^21 P23I = 2y/2C, so such state 
provides minimal S+ for a given concurrence C j 11 ' 32 The 
time tb after which the Bell inequality is no longer vi- 
olated [S+(tb) — 2] falls either into the first or second 
region, because S + (t 2 ) < 2 [it is interesting to note that 
in the third region S+ (t) passes through a minimum and 
then increases up to S+ — > 2 at t — > 00]. The time tb 
can be easily calculated if S+(t{) > 2, so that tb falls 
into the second region and therefore 

tb = (X2/2) ln(\/2 sin 2/3). (7) 

This case is realized when pure dephasing is relatively 
weak: T 1 /T 2 < ln(V2sin2/3)/[21n(4 - 2\/2)\; since 
Tx/T 2 > 1/2, it also requires sin 2/3 > 2^2- 2. [For 
Ti/T 2 = 1/2 Eq. has been obtained in Ref. [TE] 
Notice that tb in Eq. ([7]) corresponds to the condition 
C = 1/V2, while in general tb corresponds to C < l/v2 
because of the inequality-^ S+ > 2\[2C. 

Now let us focus on calculating the duration te of the 
entanglement survival, duration tb of the Bell inequal- 
ity violation, and their ratio te/tb at non-zero temper- 
ature. For simplicity we limit ourselves to the case of 
maximally entangled initial state {(3 — 7r/4), and we also 
assume equal energy relaxation, splitting and tempera- 
ture for both qubits: Tf = T-f = T lf A„ = A 6 = A, and 
9 a = 9i, = 9 (we do not need to assume equal dephas- 
ing, since it can be characterized by only one parameter 
T 2 ). As follows from Eq. the entanglement dura- 
tion te can be calculated numerically using the equation 
|p23 j = \fpi 1/344- Figure [1] shows te (normalized by T 2 ) 
as a function of the ratio T\/T 2 for several values of the 
normalized inverse temperature A/ 9. As we see, in a 
typical experimental regime^ when A/ 9 ~ 10, the ratio 
te/T 2 does not depend much on T\/T 2 when T\ is larger 
but comparable to T 2 (which is also typical experimen- 
tally). In other words, te is approximately proportional 



to T 2 , and in this regime te also has crudely inverse de- 
pendence on temperature [see Eq. ([8]) below] . 

Analytical formulas for te can be easily obtained 
in the limiting cases. In absence of pure dephasing 
(Ti/T 2 = 1/2) and low temperature (9 < A) we 
find t e /T 2 w A/26* - ln(2\/2 + 2)/2 w A/29 ~ 0.79, 
while at high temperature (9 ^> A) we have te/T 2 
ln(\/2 + l)/2 « 0.44. In the case of strong dephas- 
ing (T1/T2 ^> 1) we find (neglecting some corrections) 
t e /T 2 ~ A/(40) + ln(T x /T 2 )/2. 

However, these asymptotic formulas are not very rel- 
evant to a typical experimental situation with phase 
qubitSfi in which T\ > T 2 . As another way to approx- 
imate te we have chosen the value at the minimum of 
the curves in Fig. [TJ this minimum occurs at the ratios 
T\/T 2 somewhat close to the experimental values, and 
the result is naturally not much sensitive to T\/T 2 in 
a significantly broad range. For sufficiently small tem- 
peratures (A/ 9 > 2) we have obtained approximation 
(T E /T 2 ) min w A/46>+ln(3 3 / 4 /2) w A/4(9+0.13 and found 
that the minimum occurs at T\/T 2 k, (i£;/T2) m i n /ln3. 
So, as the crudest approximation in the experimentally- 
relevant regime {9/ A ~ 10~\ T x /T 2 > 1), the two-qubit 
entanglement lasts for (see dashed lines in Fig. [T]) 

T E T 2 A/49. (8) 

The duration tb of the Bell inequality violation is cal- 
culated using Eq. ([5]) as S+(tb) = 2. Solid and dot- 
ted lines in Fig. [2] show numerical results for tb (in 
units of T 2 ) as a function of the ratio T\/T 2 for low and 
high temperatures: A/ 8 = 15 and 0. The curves are 
almost indistinguishable, that means that tb is practi- 
cally independent of the temperature for fixed T\ and 
T 2 . Notice that each curve consists of a constant (hor- 
izontal) part and an increasing part, which correspond 
to two terms in Eq. ([5]). It can be shown that at zero 
temperature the horizontal part is realized at T\/T 2 < 
ln2/[41n(4 - 2y/2)] sa 1.1, while at high temperature 
(0 > A) it is realized at T x /T 2 < 1. The horizontal 
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FIG. 3: The ratio te/tb for the maximally entangled initial 
state and several values of the temperature 9. 

part corresponds to the first term in Eq. ((5|) dominating 
at tb'- 5+ = 2y2 exp(— 2t/T 2 ), so at sufficiently weak 
pure dephasing we have tb/T 2 = In 2/4 « 0.17 [see also 
Eq. ((?])] . In the opposite case of strong pure dephasing 
(Ti/T 2 » 1) the duration tb is the solution of the equa- 
tion t b /T 2 = ln[Ti/(4r B )]/4 (dashed line in Fig. [2D, so 
roughly t b /T 2 ~ \n{T x /T 2 )/A (dot-dashed line in Fig.©. 
Combining these results, we get a crude estimate: 

t b ^ T 2 max{0.17, 0.25 ha(T 1( /T 2 )}. (9) 

Figure [3] shows the ratio te / tb of the survival dura- 
tions of entanglement and the Bell inequality violation. 
We see that the ratio te / tb increases with the decrease 



of temperature and decrease of the pure dephasing contri- 
bution, which are both the desired experimental regimes. 
(This rule does not work in the experimentally irrelevant 
regime 9 » A and T\ < T 2 .) Notice that the kinks 
on the curves correspond to the change of the domi- 
nating term in Eq. (|5|). In absence of pure dephasing 
(Ti/T 2 = 1/2) the low-temperature result (9 < A) is 
te/tb ~ (2/ In 2) [A/0 - ln(2\/2 + 2)], while at 9 > A 
the ratio is t e /t b ~ 21n(\/2 + l)/ln2 » 2.5. In the 
limit of strong pure dephasing (T\/T 2 3> 1) the asymp- 
totic result is t e /t b ~ 2 + (A/0)/ln(Ti/T 2 ) (as we see, 
te > "2tb for any parameters). In the experimentally rel- 
evant regime when 9/ A ~ 10 _1 and T\jT 2 > 1, the ratio 
can be obtained from Eqs. §8§ and ©, giving a crude 
estimate t e /t b ~ (A/6) min{1.5, 1/ m(Ti/T 2 )}. 

For an experimental estimate let us choose parameters 
typical for best present-day experiments with supercon- 
ducting phase qubitsii A/2-kH ~ 6 GHz, 9 ~ 50 mK, 
T x ~ 450 ns, T 2 ~ 300 ns. Then A/0 ~ 6, T x /T 2 ~ 1.5, 
and we obtain r^; ~ 470 ns, ~ 60 ns, and te/tb — 7.7. 

In conclusion, we have found that in the Hilbert- 
Schmidt metric only 1.085% of entangled states violate 
the Bell inequality, thus explaining why entanglement 
can last for a significantly longer time (te) than the 
Bell inequality violation (tb)- Using the technique of 
Kraus operators, we have considered local decoherence 
due to dephasing and energy relaxation at finite temper- 
ature, and for this model calculated te, tb, and their 
ratio te /tb ■ The work was supported by NSA and DTO 
under ARO grant W911NF-04-1-0204. 
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Notice that when If / T\, the second term in Eq. (0) is 
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though the benefit is not significant if we need S+ ^2.2. 
The statement in Ref. that any mixed state with S+ = 
2v / 2C > 2 is maximally entangled, is incorrect (here max- 
imum entanglement means that C cannot be increased by 



any two-qubit unitary transformation) . As a counterexam- 
ple, consider the states p = /|*)(*| + (1 - /)|00)(00|, pro- 
duced from the initial state ((2} due to zero-temperature 
energy relaxation (T 2 = 2Ti, 6 = 0, / = e~ t/Tl ). Any two 
such states with the same / but different initial parameter 
P can obviously be connected by a unitary transformation 
(involving only the subspace spanned by |01) and 1 10} ) , 
while they have different concurrence C given by Eq. 
Finally, as follows from our analysis, there is a finite range 
of parameters / and /3, in which S+ = 2y / 2C; in this range 
the concurrence can still be varied by unitary transforma- 
tions varying (5, contradicting the statement of Ref. |ll|. 



